CHAPTER XXXVIII. 


ERRORS OR UNCERTAINTIES OF OBSERVATIONS. 


1744. Suppose a large number of observations to be made 
to ascertain the measurement of some physical element. To 
fix the ideas take one of the simplest kind, the distance 
between two marked points A and B on a straight rod. 
Suppose the distance AB to be roughly known to be 10 feet 
long, but that its true value 7 is unknown to the observers, 
of whom there are many, but known to some other person. 
And suppose that as great accuracy as possible is required. 
Out of a large number of observations by careful observers, 
it is clear that there will be none of them which differ very 
much from the true value T. The more care is taken, and 
the more accurate the means of measurement at disposal, 
the closer will the estimates be together. And it is a matter 
of experience that slight over estimates are as likely as 
under estimates, and occur with equal frequency. Absolute 
“mistakes” of counting feet or inches, or of registration of 
units, or of the use of the instruments we are not considering, 
In fact we eliminate from this explanation any errors which 
are of the class of careless “blunders.” 

It will be found by the person who knows the true value 
T, that very few of the estimates differ from T by as much 
as 4 an inch either way; fewer still by ł of an inch, still 
fewer by a whole inch, whilst errors of 4 or 5 inches would 
not occur in the tabulated results of the observations at all. 
And if the number of observations which give an error be- 
tween x and a+dz be represented graphically, it will be 
found that the graph takes the form of a curve symmetrical 
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about the y-axis, having a maximum ordinate at the origin, 
falling rapidly to the 2-axis, the ordinate speedily becoming 
insensibly small (see Fig. 586). 


bs 


(6) x 
Fig. 586. 


1745. It follows, therefore, that for the existence of an 
error of magnitude lying between x and x+dz, there will be 
a far greater probability when w is small than when z is 
large; ùe. a far greater number of errors of observation will 
fall between x and #+-dz for small values of x than for larger 
ones. Let ¢(x)dx be that number. We wish to examine the 
nature of this function g(x). And about it we know that 

(i) it decreases very rapidly as x increases ; 

(ii) it must be such as to become insensibly small within a 
short range of values of z; 

(iii) it must be an even function of g, as errors of excess or 
défect are equally numerous within corresponding 
limits ; 

(iv) it must contain some constant or constants depending 
upon the goodness of the observation, the training 
and competence of the observer, the accuracy of the 
instruments used, and the circumstances under which 
the observation is made; 


(v) the number of observations must be fe $(x)dx, and 
supposing N be this number, the chance that the 


error of any particular observation lies between 


x and «+dr=¢(x)da/N=v-(x)da, say. 


1746. Laplace’s Investigation. 
Starting with the hypothesis that an error in an observation 
is due to no one single cause, but is the aggregate of the 


www.rcin.org.pl 


THE FREQUENCY LAW. 855 


cumulative effects of a large number of causes, each pro- 
ducing its own separate effect, and that these effects are 
extremely small, and as likely to be positive as negative, 
Laplace has shown by a very laborious and difficult investi- 
gation that the chance that the error lies in magnitude 


between x and «+ dz, viz. y (x)dx, is Me e-°" dz for some value 


of w which depends upon the goodness of the observation. 
The argument is of such length that we must refer the reader 
to Laplace’s original work (Théorie Analytique des Probabilités). 
We therefore assume the law as our fundamental hypothesis 
in what follows. A good idea of the principal steps in the 
process, which avoids the obscurity of the original work of 
Laplace, will be found in Airy’s Theory of Errors of Observation, 
pages 7 to 15. Todhunter’s History of Probability, Arts. 1001 
onwards, may be consulted, also a paper by Leslie Ellis (Trans. 
Camb. Phil. Soc., viii.), and a paper by Merriman (Trans. Conn. 
Acad., iv.). 


1747. The Frequency Law. 


The law ye]? e~°* ig termed the law of “ Facility ” or 


“ Frequency” of Errors. It will be noticed at once that this 
is a probable law, for it answers all the requirements laid 
down in Art. 1745. It has a maximum at c=—0, it is an even 
function of z, it contains an arbitrary constant w, it diminishes 
with great rapidity as x increases, and speedily becomes of 
insensible magnitude, and 


ie o(a)da=N| [Semde=N. 


1748. Weight and Modulus. 

The constant w is called the weight of the observation. It 
is sometimes replaced by p Then ¢ or +. is called the 

w 

modulus. The weight w measures the care, skill and precision 
of the observer, the goodness of his instruments and the 
excellence of the conditions under which the observation is 
made. 
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1749. The ordinary method of estimating the value of a 
physical element of which a number of presumably equally good 
measurements have been made is to take the arithmetical 
mean of the result. As a matter of experience this gives good 
results, and therefore this mean is frequently adopted as giving 
the best estimate available, and regarded as the most likely 
value. If we might assume this, the above law of Facility of 
Errors easily follows. 

Let T be the true value of the measured quantity, 7 being 
unknown. Let 2,, 2,,... Zn be n independent results of obser- 
vation; (x) the law of Facility. 

Then z,—T, z,—T,...2,—T are the actual errors, some 
positive, some negative, and the à priori probability of the 
coexistence of these errors is proportional to the product 

P=¢(2,—T)¢(2,—T) ... 6(@,—T). 
Then, by the principles of inverse probability, the probability 


ihat hadran tie tean 2 wal Feat in Bae / fP dT, 


the limits being such that the integration is conducted over 
all values of T which it is capable of assuming. That is, 
after the observations were made, the probability that T is 
the true value is also proportional to the product P, and 
therefore this expression is to be made a maximum by 
variation of T. Taking logarithms and differentiating, we 


have > p (2,—T)/p(2,—T)=0. 
1 

Now, if we take for T the arithmetic mean of the observa- 
tions, this equation is to hold when nT=5 z,. To find the 
form of ¢ which will satisfy these requirements, take the case 
Zo =Z =... =n =% — nr. Then 

nT =z +(n—1) z=2,+(n—1)(z,—nr)=nz,—n(n—1)7, 

i.e. tatia odp %,—T =(2,—2,)+(4—T)=—7, 


—T=—r, ete.; 
Pa D p ntan) 
pa T) T”) pa T)? 
$(n—1) 7 $(—7) 


or 


(n—1)7¢(n—1) 7 (—7) G(—7)’ 


www.rcin.org.pl 


THE FREQUENCY LAW. 857 


which is independent of n; and this is to be true for all 
positive integral values of n. 


This will be satisfied if œ be such that  £“_const,—0; 
u $(u) 


whence log @(u)=C and ¢(u)=Ae i? 


And since ¢(u) is to decrease as u increases, C must be 


negative. Let gaii Then ¢(u)=4e “. Again, if N be 
the total number of observations, 
æ w u? 
N= | plu) du= | Ae Odu=Aelx; i A=NjeJz, 
; 2 
& = ae c* 
ie p(a)=— 64, 


which establishes the law of facility under the hypothesis 


specified as to the Arithmetic mean. 

This remark is made by Dr. Glaisher in the solutions of 
the Senate H. Problems for 1878, pages 167, 168, where there 
will also be found a concise account of the allied subject of the 
principle of “ Least Squares.” [See also Todhunter, Hist., 
Art, 1014.] 


1750. Mean of the Errors, Mean of the Squares, Error of Mean 
Square, Probable Errors. 
The following facts will now appear : 


(1) The mean of all the positive errors 


E R — G d: 
j o ON ry j c 1 
Se SS 
| otis Tř dz Jr ~ w 
0 eT 
(2) The mean of all the negative errors with their signs 
changed is also kas ea i 
pi Ja Jro 
(3) The mean of all the errors taken positively is Pero ANE ; 
alm Vw 
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(4) The mean of the re o all the errors 


aa -ag 
E =% aoe Sr he oht 
= Ro eS Pee ee TT —— 
, ae 2 2w 
—=e “dr 
Í —2 C4 T 
(5) The “ Error of Mean Square,” i.e. the square root of the 
mean of the squares of the errors, = =- This is the 


abscissa of the point of inflexion on the Probability Curve 
y=e *, 
(6) The “ Probable Error,” which is such that the number 


of positive errors which are greater than itself is equal to the 
number which are less, is given by the value of p, where 


el ATF ak, ight Msg teks 
f=- Pr inky de=5| Per oot y dx=-. 


p 
Loge 

Let z=cz. Then Fy e7” dz=0°25. 
Vr Jo 


Tables have been calculated for the values of this integral 
for various values of the upper limit [Kramp’s Refractions ; 
Encyc. Metropol., “ Theory of Probabilities ”], and interpola- 


tion from them gives A ='476948.... Hence the “ Probable 
Error ”=476948... c or 476948... /Vw. 


1751. Kramp’s Table is given by Airy (Th. of Errors, p. 22), 
also by De Morgan (Diff. Cale., p. 657). We reproduce Airy’s 
abstract of this table for convenience for other purposes. 


Integral tabulated, Z= 4 f; e” dz. 
~r Jo 


x ri x“ I £ I x I 
0°0 0°000000 1:0 0°421350 20 0497661 3°0 0499988 
0'1 0056232 ll 0°440103 2'1 0°498510 
0'2 07111351 1:2 0°455157 2'2 0°499068 
03 0164313 1:3 0°467004 2'3 0°499428 
0'4 0214196 14 0476143 2'4 0499655 
0° 0260250 15 0°483053 2°5 0°499796 
06 0°301928 16 0°488174 2'6 0°499881 
0°77 0338901 1°77 0491895 2°77 0499932 
0'8 0371051 18 0494545 2'8 0499962 
0'9 0398454 19 0496395 2'9 0499979 œ 0°500000 
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1752. Relative Magnitude of Probable Error, Mean Error, 
Error of Mean Square, Modulus. 

To sum up, we have 

Probable Error ='476948.../ Vw; 

Mean Error=1//27w='564189...//w ; 

Error of Mean Square=1//2w="707107...//w; 

Modulus=1//w ; 
in each case varying inversely as the square root of. the 
weight, t.e. directly as the modulus; and obviously, when any 
one of these is found the rest may be deduced. They are 
arranged in ascending order of magnitude. 

Taking the z-axis as the axis of magnitude of errors and 
the y-axis as the axis of frequency, Fig. 587 will exhibit to 
the eye the relative magnitude of these errors and the fall in 
frequency. The figure is that given by Airy (loc. cit. sup.). 
The abscissa is the ratio of the magnitude of an error to the 
modulus. The points P, M in the figure indicate respectively 
the abscissae for Probable and Mean Error. 

y 


to 
0-9 
os 

> 0:7 
g o6 
£05 
Doa 


Sik 


o P M M.Sq. Mod. 2 
Axis of magnitude of error x 


The Probability Curve y= Pm 
Fig. 587. 

1753. Several Observations. Resultant Weight. 

Suppose there to be a result b dependent upon two observations 
a, and a, of weights w, w, respectively, say b=¢(a,, a,). To 
find the weight of the result. 

Let xı, £, be the actual errors and z the consequent error 
in b; all being small quantities of the first order, then to that 


ð ð 
order z= Si + i L_= þa, Tı + Pato, SAY. 
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The chance of the co-existence of errors in a, and a, 
respectively between 2, and 2,+dz, for the one and a, and 
£a +da, for the other is 


@, — 2 w - 2 
qem de, . |e wt dy., 
Tv mT 


Therefore writing Ean g+ Pa and 4=> ga, the 
1 2 1 


chance of an error in b lying between z and z+ dz is 
%a, 3/2 ie | 
vasf” do E -omt -o ( 5") (5 ) | 

peek ‘eset i — e A 
T — 00 Pa, J 


a= o || Cee — Az)? 
_ Vow, e-e” dz f e opa dx 
a Pad -o 


Ves e Top Ue = ni eee Toa ee de. 
1% 


The law of facility hry the compound result ¢(&,, a,) is 
therefore of precisely the same form as that for each of the 
original peria but the weight of the combined result 

d 1 /0¢\? ; 
is w, given by i (8) Tae . And exactly in the 
same way if b depends upon several observations @,, dy, ... @n 
of weights w,, œz... œn respectively, we have a resultant 


that is, 


; r ; RE d$ 
weight w for the cumulative measure given by Ep E) ` 
It follows that, writing P.E. for Probable Error, 


[P.E. in $(a,, a, ...)P=(PE. in ap (3E 3) +(P.E. in a,)? e3) A 


and the same law of combination holds for Mean Error (M.E.) 
or Error of Mean Square (E.M.S.). 


1754. For example, if we require the weight of the Arith- 
metic Mean of n observations of equal weights w, 
* Rett ex ee 
b= 2 a/n and a WA t.e. WO=NW,. 
That is the weight of the combination is n times the weight 
of any of the original observations, and 


the Probable Error in b=(P.E. in any of the @’s)/vn, ete. 
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Similarly the weight of a resultant pa,+ga,+ra,+... is 


given by bisa ceca 
i eee A S 
w @, W2) Wg 
and if ,=w,=w,=..., (ee 
1 


1755. If observations be taken upon a single physical 
element, and the weights and probable errors of the several 
observations (a,, a, @g, ...) be respectively (w,, wz, ws, ---) and 
(ei €z» €z, ---), Whilst w and e are those of a resultant formed 
according to the law Yp,a,/Sp,, which is the usual form 
adopted, where (p,, Pz, Pz, -.-) are certain constant multipliers, 
called “combination weights,” to be so determined as to give 
a minimum probable error in that resultant, we have 


2 2 
=a (2) +en( Zt) rá 
and differentiating with regard to p,, Pa Ps -> 
f Pe =p =p =... = Epe Zr, 
Mg Pilm =p w = pw... 
ùe. the combination weights are to be proportional to the 
theoretical weights. Moreover, it follows that 
ama atat we Or w=w tw tot., 

and the theoretical weight of the result is equal to the sum of 
the theoretical weiglits of the several collateral measures (see 
Airy, Th. Err., p. 56). 

1756. To estimate the actual value of the weight of a 
series of observations upon a single physical element, we have 


seen that ee =mean of squares of the errors. 


If then the actual errors of each observation were known, 
we should have a rule to determine w. But the exact measure- 
ment of the quantity upon which the observations are made 
is rarely known, Let T be its true value, A,, A,,...A, the 
observed values. Then A,—T, A,—T, etc., are the actual errors, 
and D7 $4- But T being unknown, we have to 
approximate. Let us adopt the arithmetical mean of the 


observations as the value of T, and write T=:5, A,, which 
1 
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is known as the “apparent value,” but is not necessarily the 
true one. This gives as an approximation 


ZAP A+... + At 2TA, +4, t.)+nT =>) AP 00, 


i.e. as an approximation we have p75» a(S 4,) 
=(ean of ae I (op or of i, 
of observations of observations/ 

1757. Determination of the “ Error of Mean Square,” “ Probable 
Error,” etc., of a Measurement of an Element from the Apparent 
Errors. 

Since the true value of the measured element is rarely or 
never known, we have to devise a method of obtaining the Error 
of Mean Square, etc., by some way other than as being 1//2w, 
which would require a knowledge of w. Let A,, A,, Áz, 
be the actual results of n independent observations on the 
single physical element in question, aj, a,, a3,... the actual 
errors, T the true value; then 4,=T+a,, 4,=T+ ag, ete. 

Let M and m be the arithmetic means of the A’s and of 
the a’s. Then 


l= 1 
a,—m=Ap—T~_ y, (A,— T)=4,—7 24, =4,— M. 


The difference a,—m, viz. the difference between the actual 
error and the mean of the actual errors, is called the “ Apparent 
Error.” And the sum of the squares of the Apparent Errors 


=>) (a-—m)?=La,?— 2m .nm+ nmi=Sa,2—* (Xa,)*. 
Therefore, if @=3(4,—M)?, we have Q=Za,2—* (Za,}. 
Now let e be the error of mean square of each measure. 
Then (Art. 1750, 5) e => ay2, i.e. Sa 2=ne?. 


Again, the square of aa of error in BA, 
=(Error of mean square in }4,)}? 


=>) (Error of mean square in A,)* 
ie 
=ne (Art. 1753); 


*, sum of squares of Apparent Errors=ne—"ne=(n— l)e. 
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Hence ath toe ; and Q being known, this determines e. 


Since the Error of mean square=1/./2w, we have 


pee 1/29. 
Also Mean Error=—— = 2 Q B ; 
J= Tw mnl 

Probable Error= 0.476948 _ 0476948.. al 22... 

Jw fet 


1758. Again, since the Error of mean square of the mean 
of n independent measures of a physical quantity 


A y x Error of mean square of any one measure (Art. 1754) 


Jn 
= e= ú ET 1’ we also have 
z es 
Mean Error ja 2 Q 
of the mean}! Va n(n—1)’ 


Probable Error aye a 
of the mean \= aiei n(n—1) 


1759, Case of a System of Physical Elements. 

Suppose next that it is required to discover the values of 
a certain set of physical elements £, y, ¢,..., and that observa- 
tions upon certain connected groups of them have been taken 
giving results of the form 

plê n, &-)=Ny, pÉ m G -.-)= Ne, ete., 
the forms of ¢,, ¢,, etc., being known, and all the constants 
involved being known from theoretical or other considerations, 
whilst N,, N,,... are the results of observation, and therefore 
subject to small errors. 

Theoretically, if the number (m) of observations be the 
same as the number (u) of elements to be found, there will be a 
definite number of sets of solutions of these equations depend- 
ing upon the degrees of the several functions. If, however, 
the number of observations exceed the number of elements, 
it will not in general be possible to satisfy all the equations 
by the same values of Ê n, €, etc., and it becomes important to 
examine a method of finding their most probable values under 
the circumstances. 
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1760. Reduction of the Equations to Linear Form. 

The observed quantities N,, N,, etc., will not differ largely 
from those which would give true values to £, n, ĉ etc., and if we 
solve u of these equations we shall obtain close approximations 
to the values of £, y, & etc., or in some cases such close approxi- 
mations may be otherwise available. Let these approximate 
values be a, B, y, etc., and 2, y, 2, etc., the small residuals of 
the true values of £, n, Ẹ, etc., so that £=a+a, n=6+y, ete., 
and these residuals being small their second and higher powers 
and products may be rejected, and each equation of form 
Pilf, n §.-.)=N; may be regarded as reduced after expansion 
of (a+, 8+y, ...) by Taylor’s theorem to the type 

: azt+tby+oz+...=Ny, 
such equations being m in number. Now n; being itself the 
result of the subtraction of ¢(a, 8, y, ...) and various second 
and higher order small quantities from N, depends upon the 
observations, and is a small quantity subject to error, whilst 


Qi, Di, ¢;,... are supposed known from theoretical or other 
considerations. 


1761. The Equations of Condition. 

We therefore have m linear equations connecting u un- 
knowns v, y, 2, etc, u being <m. Let a typical equation 
be aw+by+cz+...—n=0, where i=1, 2, 3,...m. We need 
not for the moment consider g, y, Z, ... to be small. 

These m equations are not in general capable of being satisfied 
by the same values of v, y, 2,..., but we have to obtain the 
most probable values of æ, y, z,... from them; that is, as good 
an approximation as we can under the circumstances. 

These equations are called the “ Equations of Condition.” 


1762. Standardisation of the Equations. 

As to the several results of observation, n,, no, M,,.--, let 
us suppose that they are each the result of several separate 
and independent observations; e.g. taking the typical case n,, 
suppose it to have been formed as the arithmetic mean of 
w; observations upon the value of aa+by-+..., and suppose 
all these w; observations to be equally good observations. 
Then the weight of this observation is proportional to œ. 


www.rcin.org.pl 


PRINCIPLE OF LEAST SQUARES. 865 


Therefore, unless the number of observations in forming 
Ni, Na, Ng, e has been the same and the individual observa- 
tions equally good, some of the Equations of Condition will 
have greater importance than others. 

If n; be found by w; observations, each with the same pro- 
bable error e, the probable error in n; is e//w;, and the probable 
error in n.se; is e. 

Hence, if we multiply the Equations of Condition by 
Vow, Jos, Vw, etc., we get another group in which the pro- 
bable errors of the right-hand sides are each e. 

We shall suppose our m Equations of Condition to have 
been already subjected to this preparation, and therefore 
suppose that the quantities n,, na, ng, ... which occur are sub- 
ject to the same probable error e. 


1763. PRINCIPLE OF LEAST SQUARES. 
If 2, Yos %,--- be the most probable values of æ, y, z, 
respectively, then, by the nature of the case, 


Ato tbyotceot-.-— 


is a small quantity of the nature of an error. Call it v, 
Then the probability of the occurrence of the error v; being 


af ote doy, the probability of the co-existence of errors 


ae eg es 
C1, Wey heey. By 1B iee %2? dv, and as these errors have 
1 


occurred through taking x, Yo, Zo» ---, etc., as the true values 

of £, Y, 2, --., etc., the probability that £o, Yo, etc., are the true 

values is il a e°} dw; If | ae ie e dy. in 
l æ 


=o — —-o i 
which the PRUTE E. is a definite constant; and, supposing 
the Conditional Equations to have been prepared as described 
in the preceding article, the w’s occurring are all equal. 
But in any case we have to determine 2p, Yo, ete., so that 
this probability shall be as great as possible; and this will be 


m 
achieved by making >) ow? a minimum; or, if the ws are 
1 


equal, 2v?—aminimum. The method of procedure is therefore 
called the method of “ Least Squares.” 
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1764, The “Normal” Equations. 
The primary condition for a minimum is 


Dor (a; dx y+ b; dya + . va EU, 


and therefore, on equating to zero the coefficients of dx, dyo, ---, 
we have m linear equations to determine Zp, Yq, Zos «+» , VIZ 
Lw4%,=0, Ywbv,=0, Ewcia, ete; 
or in the case when the equations have been prepared before- 
hand, so that the weights are equal, 
Lawv,=0, d0,=0, Tev,=0, ete, te; 

La?. £o + Lab. yot Eac . zot ... = 2an, 

Eba . zo + EO?. Yo +Ebe . z+... = Ebn, 

Lea. zo + Ecb. yo +Uc®? . zot.. = Ten, 

etc., 

The very compact notation [ab], [aa], ete., is often used for 
Zab, La", etc., but we adopt the sigma notation as a little 
easier to write. 

These equations determine the values of 2%, Yo, etc., so as to 
give the most probable values of x, y, etc., to satisfy the 
original group of Conditional Equations in which the n’s are 
subject to small errors. 


1765. Before proceeding further, let us examine the m 
prepared equations of type aa+by+ez+...=n, from another 
point of view. 

Multiply the several equations by Pı, Po, ---, Pm and add; 
then by qı, 2, --- Ym and add; then by 7, T2, ... Tm and add; 
and so on; viz. by « groups of multipliers, m in each group. 
We obtain u equations, 

LILA p+ Yb, Pi +HZ pit ... =EN Pi 


which are known 
as the. “ Normal” 
Equations. 


Eai +YLbiqi t+ 2Ecg:i + -= Eni (1) 
Lari +ylbhr,+2rer, t= Ennii | 
etc. 
Again multiply these by ),, Ao, ..-, A, and add, and choose 
the d’s so as to remove the terms y, 2, ..., t.e. 


AZG Pit Ndegi +Agdueyr; +... ==0, 


Arb, pi + À 220.9; HA Sbr + ss. = 0, 
ete. 
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where kiA ,p;+AoitAsrit--.; whilst Ebik;=0, Eck;=0, etc., 
and the new constant multipliers k,, ka, k3, ... replace the 
p's, 7's, 7's, etc., and d’s. 

Let w be the weight of each of the observations ny, Nng ... Nm, 
by supposition prepared to be of equal weight, and let wr, wy, 
wz -.. be the weights of the deduced values of x, y, z, 


SS ee A 
ws (Lak)? w 
And if e be the error of mean square, or the probable error 


in each of the n’s, and ez, ey, €z, --. the resulting error of mean 
square, or the probable error in the deduced values of æ, y, z, ..., 


Then 


2 
we therefore have ex = SA e, and we have to make this 


error of mean square, or this probable error, as small as 
possible with the conditions 2b,k;=0, Lek,=0, ete. 


1766. To do this we have the k’s at our disposal. Their 
number is m and their connecting equations number »—1, 
which is < m. It will be observed that the expression es 
contains only the ratios of the k’s, and when their ratios to any 
particular standard k have been fixed e, becomes determinate. 
We shall therefore in no way alter the value of e, by the 
addition of some one additional linear equation amongst 
the k’s. For convenience we take that relation as La,k;=1, 
which will give ~= Xn, We then have to make e= Ek}? . ea 
minimum with the u conditions La,k,==1, Db,k;=0, Dc,k;=0, ete. 
We obtain at once Lk,dk,=0, Ladk,=0, Lb,dk,=0, etc., and by 
Lagrange’s method of undetermined multipliers 


k,=Aa,+Bb,+..., ka=Aa,+ Bb,+..., -- Kem =AGm+Bb,, + ..., 
whence Eke = Aak; =A. 


Also Aa? +BZab+0ĒZac+ ...=2a;k;=1, 
A£ba + BED? = CEbe + oon = bk; =0, 
AXca+BEcb + Cre? +...=Xek, =0, 

ete., 


Be A Ua e 
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whence A=| Xb?, Ebe, ... La’, Lab, Sac ... | and is known, 
Ecb, Be ... Dba, 20% Lhe ... 


Lea, Zeb, 20 ... 


and A=3k?. Therefore e,2=Ae and ¢,==e/A; and A is 
essentially positive, being the sum of a number of squares of 
real quantities. The weight of the deduced value for æ is 


v=o .w, and if we take w as unity, w= 5" 


1767. The symmetry of the work shows that the same 
process will give us a minimum error of mean square, or a 
minimum probable error also for y or for z, ete., and that the 
weight of y so deduced may be found by solving equations of 
the same form as those in group (3), but with the 1 now 
replaced by 0 in the first equation and the 0 by 1 in the 
second ; and so on for the weights of z, etc. 


1768. Again it will be noticed that if we choose our 
preliminary multipliers, viz. the p's, q's, 1's, etc, as the 
coefficients of the original prepared eonditional equations, viz. 
P= li T=, Ti=t, ete, we have k,=)a;+A2b;+A3¢;+..., 
and for this choice 

Lh2=Z(Ayai+Agd; +... kj =AyVa,k, + AQUA, +... =A\V=A. 

That is, substituting for the p’s, q’s, 7's, ... in equations of 
group (1), the equations which will give a value of æ with 
the least error of mean square, or least probable error for x 
are the “normal” equations arrived at in Art. 1764.otherwise, 
and the symmetry shows that the values of y, z, ete., will 
also be determined by the same equations with the least error. 
But as these equations are the same as those arrived at by 
making L(ae+by+...—n,)? a minimum by variation of 
x, Y, Z, ..., this is a convenient way of reproducing the equations 
for these unknowns. And the result is the same as that 
arrived at in Art. 1764, the weights of the several observations 
having been made equal by preparation of the conditional 
equations. 

1769. If the conditional equations are left unprepared, we 
arrive at the proper equations for the values of æ, y, z, ete., by 
making Lw,(a0+by+...—;)? a minimum. 
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1770. Reality of VA. 
The determinants occurring in Art. 1766 are essentially 
positive. For such a determinant as 


Za? Lab, Yac,...| occurs in squaring |a,, pg, ... Gm 
Eba, Eb, Ebc, ...| the rectangular array | b, b,,... b 
Dea, i2004 Denu, Ch) Ogee, 


in which the number of rows (u) is less than the number of 
columns (m), and is therefore expressible as the sum of the 
squares of all possible determinants which can be formed from 
the array by taking u columns (Burnside & Panton, Th. of Eq., 
p. 260). Such a determinant is therefore essentially positive. 


1771. To complete the theory we must examine how the 
quantity e is to be found from the details before us; that is, 
we are to do for the case of measurements upon a system of 
physical elements what was done in Art. 1757 for the measure- 
ment of a single element. We have used e indifferently in 
Art. 1765, etc., for either the error of mean square, the probable 
error or the mean error. We shall now define the letter as 
standing definitely for the “error of mean square” in the 
measure of an observation. Let v; be the residual error in 
ae+by+...—n,, when the values £o, Yo, Zo, ... obtained from 
the “normal” equations have been substituted for h TIEA 


Then we shall show that the equation cn ZE replaces 


that of Art. 1757. 
Let the true values of x, y, z, ... be y+dx, Yo tdy, %+62, 
etc., and let 
alto + dz) +b; (Yot+dy)+...—n=u,; (i=1 to i=m). 
Multiply by a; and add the system. Then 
La*.xy +Zab. yy Hac . Zo +... — Pan 


+2Za?. ôx + Zab. dy+ ae. dz+... = Dau; 

Da’. du + Lab. dy+ Lac. z+... = Dau. 

Similarly ba. dv+X0?. dy + Ube. z+... = bu, 
Eca. dx + Zeb. dy +E. dz +... =Xeu, etc., 


which, as in Arts. 1765, 1766, give da= ku. 
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1772. Equations of type at) +byo+ ...—n=v,(i=1 toi=m), 
multiplied by v; and added, give Zv? = —Xn,v,, since 
Lav=0, Ebv=0, Yev=—0, ete 
And equations of type a,(z,+6x)+b,(y)+ dy)+...—nmi=% give 
in the same way Luy;,= — Env. 
Hence 2u = Duy, = — Env. 


1773, Equations a)+by,+¢z,+...—n,=v;, multiplied by 
u; and added, give 


LAU; . Lo t+ Ebiti. Yot- — Enu; = Ev= =r. 


Equations a,(z.-+6x)+;,(y)+dy)+...—nm=u,, multiplied by 

u; and added, give 
Lal; . Ly + Xu, . Yo +- — Inu; 
+ au; 62+ Dbu; . ôy + ... = Sua, 

Hence Zu? =Zv? + Lau, .da+ Ubu;.dy+Zeu,. dz+.... 

And, since £u}? is the sum of the squares of the true errors of 
the observations, Du,2—me?. 

Now, in the terms Za;u;. cx +Eb;u;. y+ ..., we must neces- 
sarily approximate. H 

Take for them their mean values. Then 

Lau; . OL=(AyUy+Agla+Agugt ...) (ktit katot kytt...) , 

whose mean value is that of akyu, +agkyto -+ agkgug?+...; for, 
remembering that the errors u,, Us, Us,... may have either 
sign, all products involving errors with unequal suffixes will 
disappear in taking the mean. And the mean values of 
Uy", Ug”, Us’, ... are each e. 

Hence Ya,u,. ôx will be replaced by La,k,. e, that is e. 

Similarly Zb;u;. ôy, Xeu,;. dz, ... will be replaced by e. 

Therefore me? = Xv?+ pe’, u being the number of unknowns 
Si Ys By vee: 


2 
Hence f= Zot i 


m—p 


1774. If there be but one unknown, t.e. when the observa- 
tions are made upon a single physical element, we have 


Zur (Art. 1757.) 


R 
m—l' 


€ 
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1775. Effect of Exact Co-existent Relations. 
If, in addition to the m conditional equations of type 
az+by+...—n=0, 
there be p (<u) exact equations of type 
azt+Byt+...—y4=9, 
these latter equations may be regarded as determining p of the 
unknowns in terms of the other .—p. Upon substitution of 
these in the conditional equations, we have a system of m 
conditional equations amongst ~—p unknowns. Hence the 
error of mean square e will in this case be given by 
mya where v; is, as before, a,%)+by)+...—m, and 
the summation is from i=1 to =m. 
If u be large, or if there be several exact equations, a 
different process is usually employed to reduce the labour of the 
elimination. (For this see Chauvenet, Astron., p. 552, Vol. II.) 


1776. Finally, if ez, e,, €z, --- be the errors of mean square 
in Zp, Yo, %,---, and if X, Y, Z,... be the respective weights 
€ 


of Tos Yos Zo eee y then Er a, Ey = FS 
i JX JY 
X, Y, Z,... are to be determined as follows (Art. 1766): 


, ete., and the values of 


1 1 1 
For X: <a’: X +2ab- yt Zac: y +...=I1, 
A n awe 
Eba ¥ +b y tbe z t- =0, 
Zoa: +Ecb: yt ke pte =0, 
ete. ; 
For Y: a’: LS EEE, WN =0 
: y Y gt =0, 
i 2, 1 „d mi 
Sin’ yt bob: F+. pte =0, 
ete., 


the accented unknowns of each group not being required ; 
and such equations may obviously be written down from the 
normal equations. 
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Hence we obtain X, i.e. the value of sa (Art. 1766), ete. 


Moreover, in cases where the values of 2», Yo, Zo,- are 
expressed in terms of letters and not numerically, their 
weights may be obtained more readily, as in Art. 1753, by 
differentiation. 

This completes the details of the process to obtain the 
Mean Square error for each element, and the Probable and 
Mean errors may be at once deduced. 


1777. Order of Procedure. 

To sum up, the order of procedure is as follows: 

I. Given the m conditional equations amongst u unknowns 
(m> u) of type az+by+ez+...—n,=0, let each have been 
prepared by multiplication by the square root of its weight, 
viZ. /w;. 

II. From these prepared equations, or by differentiating 

L(ae+by+...—n,), 
deduce the normal equations and find a, Yos 2; =». 
III. Form Yv2f= 2 (amot bäYot ---— n). 
IV. Find e, the error of Mean Square of an observation 


Dv? 
from ej] 2, 
> M— u 


V. Then to find ez, ey, &, etc., in the normal equations 
replace Lan, Ybn, Yen, ... by 1, 0, 0, etc., and solve for x, say 
ay then replace Lan, Lbn, Len, ... by 0, 1, 0,..., etc., and 


solve for y, say y=» and so on; then X, Y, Z,...are the 


several weights of 2, Yo, %,-.., and the errors of Mean 
€ € 
JE’ a= Ta 
These values may also be obtained by Art. 1753 without 
the trouble of solving the normal equations when the results 
of the observations are given in letters instead of numerical 
quantities. i 


Square are e,= 


VI. Having found e, ez, ey, €,-..., we may then deduce the 
Probable Error or the Mean Error by Art. 1752. 
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1778. For further information, the reader may consult the appendix 
to Vol. II. of Chauvenet’s Sph. and Practical Asironomy. 

For those interested in the Bibliography of the subject, reference may 
be made to 

Legendre, Nouvelles Méthodes pour la détermination des orbites des Cométes, 
1806. 

Gauss, Theoria Motus Corporum Coelestiwm, 1809. 

Disquisitio de Elementis Ellipt. Palladis, 1811, etc. 

Bertrand, Méthode des moindres carrées, 1855. 

Encke, Ueber der Meth. d. Klein. Quad., Berlin (Astr. Year Book, 1834, etc.). 

Laplace, Théorie analytique des Probabilités. 

Poisson, Sur la probabilité des resultats moyens des observations (Con- 
naisance des Temps, 1827). 

Bessel, Astron. Nach. (357, 358, 399). 

Hansen, Do. (192, 292, etc.). 

Peirce, Astron. Journal (Camb. Mass., Vol. II). 

Liagre, Calc. des Prob., Brussels, 1852. 
And other references have been made to the works of Airy, Glaisher and 
Merriman in the course of this chapter. 


1779. ILLUSTRATIVE EXAMPLES. 

1. Suppose O a central station on a plain, and A, B, C, D four distant 
points. Let the angles AOB, BOC, COD, DOA be respectively estimated 
_ by p, q, T, 8, equally good measurements to be a, B, y, 5; and suppose that 
after all due care has been taken a+B+y+6 falls a little short of 360°, 
say by k". It is required to find the corrections to be applied to the several 
observations. 

Suppose the true values of the several angles to be a+”, B+y”, 
y+”, d+0”. 

Then «+y+2+w=k is an exact equation. The equations of condition 
are «=0, y=0, z=0, v+y+2-k=0, which cannot be satisfied simul- 
taneously. Making px?+qy?+rz?+s(x+y+z-—k)? a minimum, we have 
px=qy=rz= —3s(e«+y+2—k)=A, say. These are the Normal Equations. 

Thus ent, pur hs and A(Fri+t)=k-4; ie. Paad 

p g r Phl N s si 
kfn kiel kgl Bs} gl 
whence z=5] 25 y=7| ye see [25 w=" / 23, 
which give the probable values of x, y, z, w. 


2. Let p observations of the zenith distance of a circumpolar star be made 
at the wpper culmination, and q at the lower. It is required to find the co- 
latitude of the place. {Arry, p. 42, Errors of Observation. | 

Let a and b be the means of the two sets of observations. Then z,=a 
and z,=b are the estimated zenith distances at the two culminations. 
And we are to find the probable error in (a+b), which would be the 
true co-latitude if the means of the observations were accurate. 
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Let w be the weight of any of the original observations, all supposed 
of equal value ; w’ the weight of $(a+6). Then 
See EEN 


w 4 pw 4qw 4w pg 
Hence if « and «’ be the probable errors of an observation and of the 


deduced co-latitude, "=i Vets «e with the same formula connecting 
the errors of mean square and the mean errors. 


3. Consider a rod, whose accurate weight is h grammes, to be broken into 
three random pieces of unknown weights x, y, z grammes ; y and z are weighed 
together l times ; z and x, m times ; x and y, n times, and the means of the 
three sets of weighings are a,b and c grammes, and all the weighings are 
equally good observations so far as is known. It is required to find the most 
probable weights of the several parts and the probable error in each. 

[Matu. Terr., 1876.] 

Here x+y+z=h, (1); y+z=a, (2); 2+x=b, (3); s+y=c, (4). 

Equation (1) is exact. The others are subject to error. Let w be the 
“weight” of any one observation. The “weights” of the means are lw, 
mw, nw. The equations (2), (3), (4) may be written h-x-a=0, 
h-—y—b=0, h-—z—c=0, and we are to make 

l(h—a—a)?+m(h—y—b)?+n(h—z—c)? 
=a minimum with condition <+y+z=h. 
Thus, /(h-x—a)dx+ +=0, dx+ + =0, whence l(h-xw-—a)=...=...=A, 


ie. Bh—(z-+y-+2)-a-b-c=A(74+ 742), ins 2h-a—-b-c=A(} ++), 


i.e. x=h—a-(2h-a—b-—c) 


m 
: Wa TALE in et ee 


If w, be the “ weight” of this gp for 2, 
eh -HUEJ EJET: oy toca cote 
w, wl\da z) + wn =, (Se ene mn+ni+in 


Now À being known exactly, 2h-a—b-c is a known error, and it is 
the only ve error, and if Q "i the “weight” of this expression 
PEFP TSI 


Hi 8 507 2 
Be aoe oe (art. 1753), and (2h-a—b-—c) (Art. 1750). The 
latter equation is the PENT. one for 2. Hence 
1 Amn 1 2lmn(m+n) 


(2h -a—b—c), -(2h-a-b-c). 


o. (mn+nl+lm) 
The probable error for v, viz. - is such that 


w mn+ni+in 


w 4769... 
= otde= + one d ee —— al 
Ne f= BY, Settee 
; tH SIET 
AA p= E70.: X ataie an 


Suppose in the same example that A were not known, but that the 
several observations are (d4, dg, ... dz), (bi, bg, --- Öm), (C1, C25 «++ Cn). 
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We then have / equations of type y+z—a,=0, m of type z+x—6,=0, 
n of type x+y—c,=0. 
Then z, y, z are to be found from making 


l m n 
Uy +2-a,P + Bet a- b + 2 (e+y-—¢,) a minimum ; 
1 1 


from which (m+n) £o + nYo + mzy = Šo, E Sop 
1 1 


| Zo, Yo, Z being the values 


n l 
nato+ (at Osetio 2e, + 2a, | which give the minimum. 


l m 
MEZo + lYo + (l+ m)zo= Za, + Xb, 
We then have as an approximation 
l m n 
l 2 (Yo +20- 4p)? +È (zo+to— br) + Z (tot Yo- Ce) 
20 lim+n j 
4. A, B, C, D are four points in order on a straight line ; AB, BC, CD, 
AC, BD, AD are measured respectively a, B, y, 5, €, ¢ times with mean 
respective measurements a, b, c, d, e, f. Find the most probable vaiue of AB ; 
and if a= B=y=S=e=¢, find its probable error. (Marta. Trier., 1878.] 
Let AB=2, BC=y, CD=z; then we are to find a minimum for 
a(x —a)+ B(y—bP+y(z—cP+8(at+y—dP+e(yt2—eP + Clatyte—f/P 
The conditions are: 
a(a—a)+3(x+y—d)+ (x+y +2—f)=0, which determine 


BY =0)+ Oty d t eg +20) 4 ((at+yte—f)=0; bys 
y(z—e)+e(y+z—e) + ¢(a+y+2-f)=0, ioe 
In the case a= B=etc., these become 
3a+2y+z=atd+f, 2æx+4y+2z=b+đd+e+f, c+2y+3z=c+etf; 
whence 
w=}(2a-—b+d—e+f); y=}(-a+2b-—c+d+e); z=} (—b+2c-—d+ertf), 
i.e. x—-a=}(-2a—b+d-e+f), u+y—d=}(a+b-—c-2d+f), 
y—b=}(-a—-2b-—c+d+e), ytz2—e=}(—a+b+c-—2e+f), 
z—c=}(-b-2e-d+et+f), wr+y+z2-f=}(at+c+d+e-2f), 
and the sum of the squares of these six expressions is, say K. 


We also have 


Erig 1 
niet sitititbs 
1 


1 1 
gT t4+tl ++ 


ERa 1l 
gee 1+, 
z atari tae 


ie Wz=2w, wy,=2w, w,=20 by (Art. 1753), or they may be derived 
as in Art. 1776. 


ERRIA K E ES R T aR | 
Now NE SNE (Art. 1773); $ aerate Ec 


whence the Mean Errors, Méan Square Errors and Probable Errors of 
æ, Y, z may be at once written down. 
[See Sol. S.H. Prob., Glaisher, 1878, p. 165.] 
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PROBLEMS. 


1. In a plane triangle the angles 4, B, C are respectively 
measured m, n and p times, and the means of these measurements 
are respectively a, B and y, and a+B+y=x+«. The separate 
measurements are equally good. Show that if a+2, B+y, y+z be 
the true values of the angles, the probable values of z, y, z are 


—npe/5, —pme/d, —mne/5, where =np +pm + mn. 


2. In the plane triangle 4BC, the side b is to be determined in 
terms of a from the measured values of Band C. Find the actual 
error in the determination of b in terms of the actual errors of 
measurement of B and C, and the probable error of b in terms of the 
probable error of any measurement supposed to be the same for 
the measurement of any angle. Show that of all the directions 
in which the side b can be drawn, that gives the probable error 
of the determination of its length a minimum for which the 
angle C satisfies the equation 

ab (2a? + 3b?) (1 + cos? C) = (a* + 7a?b? + 264) cos C. 
(Maru. TRiros. ] 

3. At Pine Mount, a station in the U. S. Coast Survey, the angles 
subtended by four Paes stations 4, B, C, D were observed 
as follows : 

AB, weight 3, 65° 11’ 52-500; CD, weight 3, 87° 2’ 24-703 ; 
BC, weight 3, 66° 2415553; DA, weight 1, 141° 21’ 21-757. 

The five points are in one plane. It is required to estimate the 
corrected values of these angles. The result is that the several 
results in the seconds should be 53”:4145, 16”:4675, 25”:6175, 


24”°5005, the degrees and minutes being unaltered. 
(CHavvenet, Astron., II., p. 551.] 
4. Taking the equations 


z—-y+2z-3=0, 4x +y +4z-21=0, 
3z + 2y- 5z- 5=0, -2#+3y+32z-14=0, 
show that (1) the probable values of 2, y, z are 2'470, 3:551, 1-916 
respectively ; 
(2) the weights of x, y, z are 24:597, 13°648, 53:927 ; 
(3) the error of mean square of an observation, i.e. of 
the numbers 3, 5, 21, 14, is 0°284; 
(4) the errors of mean square of a, y, z are 0:057, 
0:077, 0:039 ; 
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(5) the probable errors of an observation and of 2, y, z 
are respectively 0:192, 0:038, 0-052, 0-026. 
(Gauss, Th. Motus; Cuauvenst, II., p. 521.] 

5. In finding the latitude of a place by observation of two 
meridian altitudes of a circumpolar star, p observations are made 
at the upper transit, g at the lower. Taking the probable error of 
each observation at the upper transit as ¢,, and at the lower as «,, 
and all astronomical and instrumental corrections to have been 
applied, show that the probable error in the determination of the 
latitude is J/pe,? + ge,2//2pq. 

6. If the altitudes of the upper and lower transits of several 
circumpolar stars be observed and H,, H,, H}, ... be the harmonic 
means of the numbers of observations at the upper and lower transits 
for the several stars, and all observations be equally trustworthy, 
with a common probable error e, supposing all astronomical and 
instrumental corrections to have been applied, show that the probable 


error in the determination of the latitude will be peH Ee 


7. At three stations P, Q, R on the same meridian, the zenith 
distances of n, stars are observed at each of the stations P, Q, R; 
n at P and Q; n, at Q and R; n, at R and P. Itis required to 
determine the amplitude of the portion PQ of the meridian. Show 
that there are four independent modes of determining that arc ; and 
on the supposition that the probable error of each observation is 
the same and =«, show how to determine the combination weights 
of the four measures. If aoe te Serie sb show that the square 


of the probable error in the result = : < : 


8. State the criterion for the selection of the combination weights 
of n independent measures of a magnitude. Determine the probable 
error of the result in terms of the probable errors of the n measures. 

In the observation of the zenith distances of stars for the amplitude 
of a meridian divided into four sections by three stations intermediate 
between the extreme stations, a stars are observed at the first, 
second, third only; b at the second, third, fourth; ¢ at the third, 
fourth, fifth; and the probable error of every observation is «. 
Show that there are only three independent modes of measuring 
the whole arc, and obtain equations for determining the combination 
weights of the three measures. In the case where a=b=c, prove 
that the square of the probable error of the result is 10¢?/3a. 

[Martu. Tre. ] 
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9. If a, b,c, ... be the actual errors in n measures of a physical 
element, the apparent error of each measure is defined as the 
difference of each measure from the mean. 

Let Q be the sum of the squares of the apparent errors. Then 
prove that (i) the Probable error of a measure, (ii) the Mean error 
of a measure, (iii) the Probable error of the Mean and (iv) the Mean 
error of the Mean are respectively 


0°674506 esi o-797885, |2 

Naw Em hy 

0674506 [2 0797885, 
n(n- 1) Nn 


(n-1) 


10. If we have any number of sets of n observations of the value 
of a physical element, all of which are à priori supposed to be equally 
good, and if the difference between any observation and the mean 
of the set of n observations to which it belongs be called the 
apparent error of that observation, then, assuming the usual law 
of frequency of errors, prove that the mean of the squares of the 


n—1 ; 
apparent errors =——— m*, where m? is the mean value of the square 
of an actual error of observation. [Smrru’s Prizzs.] 


11. A rod of known length / is broken into four portions. The 
lengths 2, y, z, w of these portions are measured respectively p, q, r, s 
times under the same circumstances and with the same care. The 
means of-these several measurements are a, 8, y, 5. Show that the 


$ ‘ l-(a+B+y+8) /1/1 1 1 
probable length of z is a+ a Aa par a ery Ae t +5). 


12. The angles of a geodetic triangle of known spherical excess 
are measured, and the probable errors of the several measurements 
are €,, €, € respectively. It is found that the sum of the three 
measurements needs a correction of 6”. Show that if a”, 8B”, y” be 
the corrections to be applied to the angles, 


a/e,? = B/e,.? = y/«,? = ble? +e? +«,”). 
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